Phase diagram of an asymmetric spin ladder 
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We investigate an asymmetric zig-zag spin ladder with different exchange integrals on both legs 
using bosonization and renormalization group. When the leg exchange integrals and frustration both 
are sufficiently small, renormalization group analysis shows that the Heisenberg critical point flows 
to an intermediate-coupling flxed point with gapless excitations and a vanishing spin velocity. When 
they are large, a spin gap opens and a dimer liquid is realized. Ifere, we find a continuous manifold 
of Hamiltonians with dimer product ground states, interpolating between the Majumdar-Ghosh and 
sawtooth spin-chain model. 
PACS numbers: 75.10Jm, 75.10b, 71.10.Hf 



The interplay of geometric frustration and quantum 
fluctuations, and eventually broken translational invari- 
ance, in low-dimensional spin systems gives rise to novel 
magnetic phases. Examples include: spin chains with 
interactions beyond nearest neighbors, spin ladders, tri- 
angular and Kagome systems Spin ladders in par- 
ticular, have attracted much interest and several model 
problems are well understood [||, and highlight the role 
played by frustration. Spin-isotropic two-leg ladders 
with railroad geometry quite generally lead to a sin- 
glet ground state separated by a finite excitation gap 
from the first triplet states |^ . Zig-zag ladders are more 
strongly frustrated and, depending on the ratio of the 
leg to rung exchange integrals, may have gapless spin 
liquid ground states or gapped dimer states ||^. Com- 
pounds hke Cu2{C^Hi2N2)2Cl4, and CuGeOs have been 
suggested to be described as railroad and zig-zag ladders, 
respectively 

Little work has been done on asymmetric spin ladders 
where the exchange integrals on both legs differ. Only 
the extreme case where one leg of a zig-zag ladder is 
missing entirely (sawtooth or A-chain) has been solved 
|^,||. The ground state of this model is a product of 
nearest-neighbor singlets and there is a spin gap, as in 
the Majumdar-Ghosh model (t). The properties of the 
excitations in both models are different, however. Saw- 
tooth chains may describe the material YCUO2.5 [HiH- 
Here, we perform a systematic study of asymmetric zig- 
zag spin ladders. Important questions concern possible 
new phases generated by the leg-asymmetry, their exci- 
tation spectra and their physical properties as well quan- 
titative changes brought about by leg-asymmetry to the 
excitations in the more usual phases found on zig-zag 
ladders (spin Hquid, Neel and dimer states). 

We consider a Heisenberg model on the structure 
shown in Figure |l|, and represent this as a chain with 
an alternating next-nearest-neighbor (NNN) exchange 

H = j2 {Ji^i ■ Si+i + [J2 + i-iys] Si ■ S1+2} , (1) 
I 

where Ji = 1 and J2 ± 5 are the nearest-neighbor (NN) 
and alternating NNN coupling constants, respectively. 
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FIG. 1. The asymmetric zig-zag spin ladder. 

J = is the ordinary zig-zag ladder or frustrated spin 
chain and fixing in addition J2 = \Ji gives the ex- 
actly solved Majumdar-Ghosh (M-G) model 0. The 
sawtooth chain has 5 = J2 = \ Ji [||J|]- 

Following the general approach to map the quantum 
spin model to a continuum field theory Q and using the 
standard dictionary of Abelian bosonization , we ob- 
tain an effective boson Hamiltonian H = Hq + Hi with 



H, = I dx- 



Hi= dx 
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^{x) is a bosonic phase field and n(a;) its canonically con- 
jugate momentum, u and K are the effective spin velocity 
and coupling constants, including the effects of marginal 
interactions, a is a short-distance cutoff, (73 oc 1 — J2/ J2c 
is the Umklapp-scattering amplitude and gi (x S is the 
amplitude of the alternating NNN field. In this paper, 
we only consider the case Ji/2 > J2 > S, while, for 
J2 > O.5J1, quite different field theory treatments are 
required 0,0. 

For 5 = 0, the model (|l|) is well understood 
Its elementary excitations arc spinons which are gap- 
less for J2 < J2C — 0.2412 0] when frustration is ir- 
relevant and the ground state is unique, or gapped for 
J2 > J2C when frustration is relevant and the ground 
state is doubly degenerate. In the field theory (g), (||), 
this translates into the g^ term being either marginally 
irrelevant, leading to the weak-coupling Heisenberg fixed 
point (X^-^) = 1/2, (^3^^ =0) 0, or relevant, yielding a 
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strong-coupling dimer state. When the SU(2) spin sym- 
metry is broken, a Neel state or an easy-plane spin liquid 
may form. Introducing dimerization, i.e. an alternat- 
ing Ji, produces an effective confining potential between 
the spinous. The elementary excitations become a spin 
triplet and a spin singlet above the unique ground state. 
In the language of field theory, the external dimeriza- 
tion corresponds to a relevant term sin 2$ which lifts the 
double degeneracy of the dimer ground state . 

Qualitative results on the influence of the new inter- 
action (gi) can be obtained from physical considerations 
alone. Firstly, the scaling dimensions of the Umklapp 
and nearest neighbor alternation terms and 171 are 



d u 
IT 



(8) 



rfg., = 4X , dg^= K +1 



(4) 



At the Heisenberg fixed point, 173 is marginal with dg^ = 
2, while the gi-term with dg-^ = 3/2 is relevant. We con- 
clude that gi destabilizes the isotropic Heisenberg fixed 
point and the spin liquid ground state. On the other 
hand, there is no standard strong coupling theory for 
the gi-term. Usually (e.g. 33 ±00), the boson field 
^{x) locks into a constant value with small fluctuations, 
and an associated excitation gap. Such a phase lock- 
ing, however, is forbidden by the $-prefactor to the 
cos(2$)-term in Hi. Secondly, the gi-term, induced by 
the alternating NNN interaction, does not confine the 
spinous and plays a role very different from the external 
dimerization, as can be checked in the dimer state by 
comparing the sawtooth chain |^,^ with the M-G model 
. Moreover, from differences in the size of the spin gaps 
in these two models, it follows that the gi-term quite gen- 
erally counteracts the Umklapp term while an external 
NN dimerization would cooperate. 

That gi opens no spin gap despite being a relevant 
perturbation of the Heisenberg fixed point, is also cor- 
roborated by the absence of a magnetization plateau in 
our model in small mag netic fields |T|. A necessary con- 
dition for the formation of a magnetization plateau in 
the absence of modulated external fields is an alternat- 
ing component of the exchange integrals . For an 
alternating NN exchange, a magnetization plateau is ob- 
served in small magnetic fields, but alternating NNN ex- 
change, it is only observed in high fields p3 I5|. 



We now perform a perturbative renormalization group 
(RG) analysis by mapping the model on a modified (by 
gi) classical 2D XY-model 

Introducing the re- 
duced variables 1/3 — g^/iru and yi — (7i/V27ru, we ob- 
tain the linearized RG equations 



d K 
"df 
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d yi 

dl 



~yjK' + ylK^ 

(2 - 4if )j/3 + if'yi' 

(1 - K)yi - 4if2y^y3 
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under a change of length scale a ae . These equa- 
tions, and our solutions to be discussed below, can also 
accommodate anisotropy in the NN and NNN exchange 
integrals. The RG equation for the spin velocity u is 
a consequence of the anisotropy of the gi-interaction in 
the classical 2D XY-model, i.e. its non-retarded but non- 
local character in the quantum field theory (^), and has 
been discussed before in the ID electron-phonon problem 
]l8[ . For yi — 0, one finds the three generic phases dis- 
cussed earlier (i) a weak-coupling spin liquid phase 
(53 = 0, K* > 1/2) terminating at the isotropic Heisen- 
berg fixed point isT^^) = 1/2; (ii) a strong coupling 
Neel phase g^ 00; (iii) a strong coupling dimer phase 
93 —00. Moreover, u is not renormalized for yi — 0. 
Spin-rotation invariant models scale along the separatrix 
between Neel and spin liquid phases, or along its contin- 
uation into the dimer regime. 

Taking j/i, the most relevant perturbation, finite, 
changes this simple picture. Quite generally, the cor- 
rections to the RG flow of K and g^ are positive, while 
those from 53 are negative. In the spin liquid phase, 
close to the isotropic Heisenberg fixed point, the alter- 
nating NNN exchange will increase both K and ^3, and 
therefore reverse the direction of the RG flow, compared 
to the j/i = situation. The sign of this effect agrees 
with the spin gap reduction in the dimer phase when 
going from the M-G to the sawtooth model. The mag- 
nitude of the effect is a direct consequence of the scaling 
dimension dg-^^ . Figure |^ shows a family of solutions of the 
RG equations, projected on the y^ — K -plane, which have 
been linearized for the purpose of the figure to accurately 
control spin-rotation invariance (cf. below). 




FIG. 2. The scaling trajectories for yi{l = 0) = 0.001 pro- 
jected on the j/3 — SK plane. 5K = K — 1/2, and the dot 
locates the new intermediate coupling fixed point. The Neel 
state is realized in the upper left, the dimer state in the lower 
left, and the spin liquid in the right part of the figure. The star 
locates the boundary between flows to the new flxed point, 
and into the dimer regime. 
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The flow reversal is seen in much of the first quadrant. 
When K is increased sufficiently by yi , this perturbation 
becomes irrelevant, however, and the RG flow bends back 
to a weak coupling fixed hue {K* > 1, 7/3 = 0, = 0), 
describing a spin liquid with an increased easy-plane 
anisotropy. Its spin velocity u is decreased but remains 
finite. 

More interesting is the spin-rotation invariant line 
(2/3 = 4_ftr — 2 in the linearized RG). While the ini- 
tial flow towards the Heiscnberg fixed point also reverses 
its direction under the influence of yi, it reaches a new 
finite-couphng fixed point {K* w 0.578, j/3 w 0.315, yj" ~ 
±0.544), resp. {K* = 0.6, = OA,y^ = ±0.8) in the 
linearized RG. This fixed point is attractive for spin- 
rotation invariant systems, and repulsive otherwise with 
a fiow into the spin liquid or Neel phases for easy plane 
(easy axis) anisotropy. At the fixed point, the renor- 
malized spin velocity vanishes, u* = 0. This conclu- 
sion follows from Eq. (||) and is independent of the ex- 
act value of the fixed point as long as it is located on 
the RG separatrix with 1/2 < K < 00. Consequently, 
conformal invariance is broken, and we do not expect 
power-law decay of the correlation functions. Thermo- 
dynamic properties will be different from a Heisenberg 
chain: the specific heat will show nonlinear tempera- 
ture dependence, and the magnetic susceptibility will de- 
pend on temperature. This state is significantly different 
from the Luttinger-type spin fluid state of the Heisenberg 
model. The fixed-point Hamiltonian can be rewritten as 
H* = u*H{K* , g*, g^) where H is independent of u*. 
It becomes trivial, H* = 0, at the fixed point because 
u* = 0. We interpret this as our spins effectively decou- 
pling at the lowest energy scales, i.e. a kind of asymptotic 
freedom in this spin-rotation invariant ladder. 

From standard arguments, we find that the elementary 
excitations in the fixed point, spinon and antispinon, are 
still gapless. One can also provide a variational argu- 
ment, based on an effective Hamiltonian Hi, Eq. (|^). 
Minimizing its energy directly produces pairwise kink 
and antikink solutions corresponding to the elementary 
excitations with energies proportional to 33 resp. gi. In- 
cluding the quantum fluctuations from Hq is expected 
to delocalize these kinks and reduce their energies fur- 
ther. Using now our fixed point properties 51^3 = y\^u* , 
M* — !■ makes the excitation gap vanish at the fixed 
point. Also, the numerical results of Wiessner et al., [p^ , 
indicate a paramagnetic susceptibility. 

A natural question arises here: Is this intermediate 
fixed point stable against higher order perturbations? 
Our answer is positive although we do not attach any 
particular significance to the numerical values of the cou- 
pling constants, which may change as higher-order cor- 
rections are included. Its location on the RG separatrix 
between the Neel state and the easy-plane spin liquid 
is protected by spin-rotation invariance and it can be 
pushed neither back to the Heisenberg fixed point, nor 



to strong coupling, K 00, by higher-order operators. 
The first option is inconsistent with the scaling dimen- 
sions, the second one would correspond, in the fermion 
language, to long-range pairing order which is excluded 
in ID. Scaling into the dimer regime is inconsistent with 
the absence of a magnetization plateau Q . 

Care should be taken when comparing these predic- 
tions to results, e.g., from exact numerical diagonaliza- 
tion on small clusters (size N) . The finite system size will 
stop the RG fiow at a scale In = In N, and the fixed point 
(l 00) is not reached. However, Eqs. (^)-(||) predict 
a spin velocity decreased significantly by the alternating 
component of the NNN exchange, with an unusual size 
dependence which can be evaluated by integrating the 
RG equations up to In- Due to the scale dependence 
of the renormalization, the finite-size spectrum likely ex- 
hibits significant nonlinear corrections which, again, are 
size-dependent. The velocity renormalization also sug- 
gests a scale-dependence of the magnetic susceptibility. 
A direct solution of our model in a magnetic field is rather 
involved, however, because higher-order terms lead to a 
field-induced generation of new, relevant operators. This 
will be reported later. 

Of course, when J2 increases beyond a critical value 
J2c{S) now depending on the exchange alternation, the 
RG flows to a strong coupling fixed point, which corre- 
sponds to the quantum dimer phase. For yi = 0.001, this 
critical point is indicated in {K, ?/3)-coordinates in Figure 
2 by a star. For small 6 and J2 > J2c{S), our RG equa- 
tions show that the spin gap is decreased by increasing 
d, but the system will remain in the universality class of 
the dimer liquid. 

However, we do not expect that the field theory de- 
scription is very precise for J2 near the M-G point 
J2 = O.5J1. The main reason is that the correlation 
length, which is proportional to the inverse of the energy 
gap, decreases quickly when J2 increases much beyond 
J2c{S)- To access this limit, we now discuss the influence 
of NNN exchange alternation on the ground state and 
excitations of the Majumdar- Ghosh model. For the M-G 
model 0, the two linearly independent ground states, 
say, left or right dimer ground state, are products of 
nearest-neighborly singlets, respectively 

i<i'L>= n 1*^)- n [M+i], (9) 

I — odd I — even 

where = - [i]i[]]j)/V2 denotes the sin- 

glet combinations of spin i and j. (^) also represents the 
degenerated ground state of the sawtooth model 
How are these models connected? We notice that NNN 
exchange alternation does not modify the product states 
of nearest-neighbor singlets 

Hs I ^lm) = ^(-l)'<5S,S,+2 I ^l,r) - 0. (10) 

Furthermore, we can prove that introducing NNN ex- 
change alternation into the M-G model, or changing it 
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in the sawtooth chain, will not affect their ground states. 
There is thus an entire manifold of Hamiltonians, pa- 
rameterized by S with Ji — 2J2 fixed, with doubly de- 
generate NN-dimer product ground states |$l,_r). This 
kind of ground state attracted much attention recently 
for the experimental realization of SrCM2 (-803)2 as 2D 
Shastry-Sutherland model [^9|j20|. 

While in these models, the focus is on the ground state, 
we also can characterize the elementary excitations as 
kinks and antikinks with finite excitation energies and 
different dispersions, in the dimer state of our ladder. 
Starting from the M-G model, the kinks can be thought 
of as domain walls separating the different dimer ground 
state configurations. From symmetry considerations, the 
kink and antikink properties are identical in the M-G 
model. With alternating NNN interaction, the symme- 
try between legs is broken, some properties will differ be- 
tween kinks and antikinks, in particular the dispersions. 
However, they still survive as the elementary excitations 
of the asymmetric spin ladder system. The alternating 
NNN interaction does not serve as a confining potential 
for the kink and antikink, but changes the energy gap of 
the excitations. Based on the cluster variation method 
1^,^, our numerical estimates indicate that the gap de- 
creases from 0.234 in the M-G model to 0.219 in the 
sawtooth chain with increasing S, while the ground state 
remains invariant. Particularly, for 6 = J2, i.e. the saw- 
tooth chain, the kink excitation is exactly a single spin 
on odd site and dispersionless, while an antikink is still 
a domain wall propagating with an effective mass |^,^ . 

The preceding results can also be understood from the 
point of view of the corresponding field model. The pres- 
ence of gi term 92:$cos2$ does not lift the degeneracy 
of the dimer ground state, which is the constant solution 
minimizing the energy of the Hamiltonian with 93 term. 
Near the strong coupling dimer fixed point, the term 
is much more relevant than the gi term. In this case, the 
soliton solutions of the cos 4$ sine-Gordon equation will 
survive. However, the gi term, although it is much less 
relevant than the (73 term, gives different masses to the 
sine-Gordon kink and antikink solutions. This give us a 
rough explanation to the asymmetric kink and antikink 
excitations in the presence of alternating NNN coupling. 
The difference in (numerically determined) gap size be- 
tween the M-G and sawtooth chains, is consistent with 
our RG results on the influence of 6 for smaller J2, on 
the spin gap magnitude. 

Spin-isotropic, asymmetric zig-zag ladders are de- 
scribed, in the limit of weak frustration, by an 
intermediate-coupling fixed point with gapless excita- 
tions and a vanishing spin velocity, likely indicating a 
decoupling of the spins at low energy scales. For larger 



frustration, a more usual dimer liquid phase is realized 
whose spin gap decreases with increasing leg-asymmetry. 
A continuous manifold of Hamiltonians with the same 
singlet product ground state interpolates between the 
Majumdar-Ghosh model and the sawtooth spin chain. 
In addition, we find gapless spin liquid and gapped Neel 
states with easy-plane and easy-axis anisotropy. Exten- 
sions we currently consider include external fields and 
doping with charge carriers. 
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